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^] | Abstract. We present the analysis of the first observations of the Cygnus region by the SPI spectrometer onboard 

the Integral Gamma Ray Observatory, encompassing ~ 600 ks of data. Three sources namely Cyg X-l, Cyg X-3 
and EXO 2030+375 were clearly detected. Our data illustrate the temporal variability of Cyg X-l in the energy 
£f} , range from 20 keV to 300 keV. The spectral analysis shows a remarkable stability of the Cyg X-l spectra when 

C^*) ■ averaged over one day timescale. 

The other goal of these observations is SPI inflight calibration and performance verification. The latest objective 
has been achieved as demonstrated by the results presented in this paper. 

Key words, gamma rays: observations — X-ray stars: individual: Cyg X-l, Cyg X-3, EXO 2030+375 — Black hole: 



physics — Space telescope: INTEGRAL: SPI 

1. Introduction of a blue supergiant (HDE 226868) and a compact com- 

panion generally believed to be the prototype of a black 

The spectrometer for INTEGRAL (SPI) is one of the two hok candidate . Its high energy emission is characterised 

rrj ■ main instruments onboard ESA's INTEGRAL observa- u • u . r , , . i £ -it a 

uu , by a variability on time scale ranging from milliseconds 

tory launched from Bakonour, Kasakhstan, on October tQ months Cyg x _ 3 ig & bright high _ mass x _ ray binary 

17, 2002. The aim of SPI is to perform the high resolution ^ & 4 g hf Qrbital lying &t & distance Qf _ 1Q 

spectroscopy as well as imaging of astrophysical sources in kpc (Dickey lg83) Thfi exact nature (bl&ck hole Qr neu _ 

the energy range between 20 keV and 8 MeV The imager trQn ^ q{ thg comp&ct object ig ^ ^ question Xhe 

IBIS (15 keV-10 MeV) (Ubertini et al., 1997) is comple- companion is most likely an early W olf-Rayet star (Van 

mentary to SPI with its high angular resolution. JEM-X, Kerkwijk et ^ 1992) The source is highly yariable from 

an X-ray monitor (Westergaard et al., 1997), and OMC, radio tQ hard x _ rayg ft ig Qne of thfi dozen GaJactic bina . 

an optical monitor camera (Gimenez et al, 1997), com- rieg considered as a " microquasar " (Distefano et al, 2002). 
plete the INTEGRAL observatory. 

The INTEGRAL observations of the Cygnus region EXO 2030+375 is a transient accreting X-ray pulsar 

started on 2002, November, 15 and lasted up to 2002, discovered by EXOSAT during a giant outburst in 1985 

December, 24. Three high energy sources, namely Cyg X- (Parmar et al, 1989). Localised at - 3.6° from Cyg X-3, 

1, Cyg X-3 and EXO 2030+375 were clearly detected by i* is a hi S h mass binary system consisting of a neutron 

gpj star and a Be companion. Its X-ray emission is charac- 

Cyg X-l is one of the brightest and most extensively terised h y the 41 - 8 s neutron star pulsation coupled with 
studied gamma-ray sources in the sky. The system consists recurrent outbursts occurring every 42 days and lasting 7 
to 9 days, thought to mark the periastron passage of the 
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by BATSE (Laycock et al., 2003) while a multiwavelength 
study was performed by Wilson et al. (2002). 

2. Instrument 

The SPI spectrometer consists of an array of 19 actively 
cooled high resolution Ge detectors with an area of 508 
cm 2 and a thickness of 7 cm. It is surrounded by a 5 cm 
thick BGO shield. The detectors cover an energy range 
of 20 keV - 8 MeV with an energy resolution (FWHM) 
ranging from 2 to 8 keV as a funtion of energy. In addi- 
tion to its spectroscopic capability, SPI can image the sky 
with a spatial resolution of about 2.6° over a field of view 
of 30°, thanks to its coded mask. Despite such a modest 
angular resolution, it is possible to locate sources within 
few arcminutes. See Vedrenne et al. (this issue) for more 
details. 

A complete INTEGRAL orbit lasts ~ 3 days, but sci- 
entific data cannot be accumulated when the instrument 
is crossing the radiation belts, reducing the useful time 
by 30%. A detailed description of the instrument in-flight 
performances is given in Roques et al. (this issue). 

3. Observations 

Each orbit consists in several exposures (typically, corre- 
sponding to the observation of a given source). Owing to 
the small number of detectors, imaging relies on the pos- 
sibility of performing the observations in dithering mode 
(Vedrenne et al., this issue). The pointing direction varies 
in step of 2° with a 5 x 5 quadratic or 7 hexagonal pattern 
within an exposure. Each pointing has a duration gener- 
ally fixed to ~ 30 minutes. 

Observations can be performed either in spectral 
mode, where the single events are accumulated in one 
unique spectrum for each pointing, or in photon-photon 
mode, where all events are time tagged with a temporal 
resolution of 102.4/is, and stored individually. 

The data used in this analysis were recorded from 2002 
November 15 to 2002, December 24 (table 1). Variations 
of the temperature of the Ge detectors produce a small 
gain drift in their associated electronics. Such an effect is 
easy to correct revolution per revolution, through the en- 
ergy calibration of each detector using known background 
induced lines. For the analysis described in this paper, we 
used only the single detector events (see Diehl et al., this 
issue). 

4. Data analysis 

4.1. Images analysis 

4.1.1. Misalignment between SPI and the star trackers 

The first images of both Cyg X-l and the Crab nebula 
yielded positions on average 8 arcmin away from their cor- 
rect values. Such a misalignment between the star trackers 
and the SPI axes had been measured on ground and later 



found stable after launch. A correction has thus been com- 
puted and included in the data analysis. This correction 
is similar to that obtained by Walter et al. (this issue), 
which was not available at the time of this analysis. 

4.1.2. Background determination 

The background substraction method depends on the 
pointing mode. For exposures in dithering mode, when a 
bright source illuminates a different combination of detec- 
tors in each step of the dithering pattern, it is possible to 
determine the background level for each detector indepen- 
dently. In the staring mode, the unique pointing direction 
makes it impossible to determine background mathemat- 
ically. A model of the background distributions on the 
detector plane has to be assumed. Analyses of "Empty- 
field" observations have shown that, for a given energy, 
the relative count rates of the 19 Ge detectors are broadly 
constant while the global amplitude (normalisation fac- 
tor) varies with time. We thus fix the relative count-rate 
distribution (depending on the energy band) leaving the 
normalisation as the only free parameter. 

To limit the number of free parameters, we fix the 
background values over a timescale as long as possi- 
ble. Considering the flux variations measured by the 
Anticoincidence Shield or Ge saturated events, we assume 
that the background parameters remain constant over typ- 
ically half an orbit, in both staring and dithering modes. 

4.1.3. Images 

We have used an iterative source removal algorithm to 
build sky images (SPIROS, Skinner et al., this issue). It 
is always difficult to build images with a strong and vari- 
able source. The main uncertainties come from the back- 
ground determination. Although the background is gen- 
erally quite stable, (solar flares or radiation belts regions 
excepted) , variations of several percents of the background 
along with source variability can yield unreliable results. 
In order to limit these effects, we have not used pointings 
for which reconstructed images (convolved by the reponse 
matrix), when compared to detector counts, give an un- 
acceptable x 2 ' Those precautions reduced the number of 
pointings used by ~ 10 %. For each image, we are look- 
ing for a maximum of 10 sources or all sources above a 
threshold level of 2 cr. 

Fig. 1 shows an image of the Cygnus region, in the 
energy range between 20 and 84 keV, obtained by SPI 
using all the cleaned data, i.e. ~ 600 ks. 

It demonstrates, first of all, that our misalignment cor- 
rection is valid: the SPI position for Cyg X-l turned out to 
be a 20 oo = 299.559±0.016 £ 20 oo = 35. 202±0.012, compat- 
ible, within less than 2 arcmin, with its known position. 
The statistical significance of the detection of Cyg X-l is 
- 680 cr (20 - 84 keV). 

Fig. 1 shows also the presence of Cyg X-3 at a2ooo= 
308.222 ± 0.090 £ 20 oo=40.929 ± 0.062, compatible with 
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DEC (degrees) 



Revolution Pointing 
number mode 



Date start - end 
Year /month/day - hh:mm 



Accumulation 
mode 



Useful duration 
kilosecond 



11 


Staring 


02/11/15 - 


07:17 - 


02/11/18 - 


07:02 


photon - photon 


27.8 


16 


Staring 


02/11/30 - 


06:08 - 


02/12/03 - 


05:54 


onboard spectra 


56.9 


17 


Staring 


02/12/03 - 


05:54 - 


02/12/06 - 


05:39 


photon-photon 


41.5 


18 


Staring 


02/12/06 - 


05:39 - 


02/12/09 - 


05:24 


onboard spectra 


47.4 


19 


Dithering 


02/12/09 - 


05:24 - 


02/12/12 - 


05:11 


onboard spectra 


75.7 


20 


Dithering 


02/12/12 - 


05:11 - 


02/12/15 - 


05:00 


onboard spectra 


48.5 


21 


Dithering 


02/12/15 - 


05:00 - 


02/12/18 - 


04:48 


onboard spectra 


58.0 




Dithering 


02/12/15 - 


05:00 - 


02/12/18 - 


04:48 


photon-photon 


39.3 


22 


Dithering 


02/12/18 - 


04:48 - 


02/12/21 - 


04:34 


onboard spectra 


123.2 


23 


Dithering 


02/12/21 - 


04:34 - 


02/12/24 - 


04:18 


onboard spectra 


75.8 
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Fig. 1. An Image of the Cygnus region using all data listed Fig. 2. An image of the Cygnus region as above between 
in table 1 between 20 and 84 keV. The significance levels 84 and 308 keV. The significance level for Cyg X-l is 247 
for Cyg X-l and Cyg X-3 are 684 a and 74 a respectively. G while Cyg X-3 is not detected. 



the source theoretical position, with a significance of 74 
<j. 

When selecting higher energy ranges, as in figure 2 
and 3, reconstructed using photons in the 84 - 308 keV 
and 308 - 615 keV energy ranges, only Cyg X-l continues 
to be detected. 

Figure 4 contains only the period of activity of EXO 
2030+375 , namely revolution 19 to 22, with a total of 345 
ks of observing time (see section 4.2). 

EXO 2030+375 is found at a 20 oo=308.02 ± 0.10 
^2000=37.55 ± 0.08, compatible with its known position. 



In figure 4 (20 - 50 keV), the statistical significance of 
sources is as follows: Cyg X-l is detected at 570 a, Cyg 
X-3 at 86 a and EXO 2030+375 at 36 a. 

All the other excesses are found with a significance 
less than 15 a, i.e less than 3% of the strongest source. 
The presence of noise at this level can be explained in 
terms of Cyg X-l variability, background determination 
procedure, and/or systematic effects due to uncertainties 
in the calculated coded mask response. 
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Fig. 3. An Image of the Cygnus region using all data listed 
in table 1 between 308 and 615 keV. The significance level 
for Cyg X-l is 16 a. 
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RA (degrees) 



100 200 
Significance 



Fig. 4. Image of Cygnus region for revolutions 19 to 22 
(see table 1) in the energy range between 20 and 50 keV. 

2C ^8 KeV 'ght cu-ve 



4.2. Light curves 

The light curves of the three sources are obtained simul- 
taneously, with the background fixed to the mean orbit 
values (see section 4.1.2). The source positions are left 
free but they are required to be constant for the whole set 
of observations. 

For Cyg X-3 and EXO 2030+375, the pointings have 
been regrouped to get a timescale of ~ 5-6 hours. 

The light curve of Cyg X-l (fig 5) shows a marked 
variability at low energy (20 - 48 keV), by a factor of ~ 2. 

Fig. 6 and 7 show the light curves of Cyg X-3 and EXO 
2030+375 between 20 and 45 keV. Although the statistics 
are much less impressive than in the case of Cyg X-l, we 
can see that Cyg X-3 is detected through all the revo- 
lutions, while EXO 2030+375 exhibits temporal evolution 
corresponding to an expected periastron outburst, with no 
detection before day 2002/12/09, then a significant emis- 
sion from 2002/12/10 to 02/12/20 (middle of revolution 
19 through revolution 22, see table 1). 

4.3. Spectra 

Counts and photon spectra for all the sources in the 
field of view are extracted simultaneously by an algorithm 
which fits all the positions and intensities of a set of given 
sources. To build photons spectra, we need the energy re- 
sponse matrix corresponding to each pointing. We used 
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Fig. 5. The light curve of Cygnus X-l between 20 and 48 
keV and 48 and 142 keV with a time resolution of ~ 30 
minutes 

2003, May version of the energy response matrix (Sturner 
et al., this issue). 

4.3.1. Cyg X-l 

Fig. 8 shows the spectrum of Cyg X-l averaged over the 
total exposure time of ~ 600 ks. 

The source is detected up to ~ 1 MeV and its spec- 
trum is well described, at least up to ~ 300 keV, by the 
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45 keV light curve 
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Fig. 6. Light curves of Cygnus X-3 with a time resolution 
of 5-6 hours (20-45 keV) 

20 - 45 keV light curve 



40 50 60 

Day from the 17th of October, 2002 



Fig. 7. Light curves of EXO 2030+375 with a time reso- 
lution of 5-6 hours (20-45 keV) 



usual Sunyaev-Titarchuck (S-T) inverse Compton emis- 
sion model (Sunyaev and Titarchuck, 1980) or by a power 
law with an exponential cutoff. For a spherical geometry, 
the temperature (kT e ) is ~ 48 keV and the optical depth 
is r ~ 2.5 for the S-T model while for the second model the 
power law index is ~ 1.5 and the cutoff energy ~ 155 keV. 
A thermal law gives a temperature of ~ 70 keV. In oder to 
compare our results with those obtained by OSSE/GRO, 
where the fitting was performed between 60 keV and 4 
MeV (Phlips et al., 1996), we computed also the best fit- 
ting parameters above 60 keV. The new fit yielded higher 
kT e ~ 58 keV with a lower optical depth r ~ 2.0, a set of 
parameters fully compatible with those of OSSE/CGRO. 

Above 300 keV, we detect an excess with respect to 
the S-T model, similar to that measured by OSSE/GRO 
(Phlips et al., 1996). We cannot exclude that the variabil- 
ity of the source produces such a feature when the data 
are integrated over a long time. 
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Fig. 8. Mean photon spectra of Cyg X-l for the entire 
set of observations (~ 600 ks). Solid line is Sunyaev- 
Titarchuck (S-T) model and dotted line is the power law 
+ exponential cutoff model. Upper limits are at 1 a. 

The individual spectra, built on timescale of one or- 
bit, reflect the variability of the source intensity, but the 
spectral parameters remain remarkably stable, with kT e 
varying from 45 to 54 keV and r from 2.2 to 2.6, the power 
law with an exponential cutoff model gives a between 1.5 
and 1.7 and a cutoff energy in the energy range 144 - 205 
keV. We have to note that the x 2 ? obviously rather high, 
are lower with the exponential cutoff power law model. 

4.3.2. Cyg X-3 

The Cyg X-3 average spectrum is shown on fig. 9. 



10 
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100 
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Fig. 9. The average spectrum of Cyg X-3 for a total ex- 
position of 600 ks. The solid line is a best fit power law 
model, with an index of 3.45 ± 0.01 and a flux at 50 keV of 
2.19 ± 0.01 10" 4 photons-cm- 2 -s- 1 -keV _1 . Upper limits 
are at 1 a. 



Cyg X-3 is detected up to ~100 keV and its av- 
erage photon spectrum is well fitted by a power law 
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of index -3.45 with a flux at 50 keV of -2.2 x 10 -4 
photons-cm -2 •s _1 -keV _1 . This index value is similar to 
that reported by OSSE (McCollough et al., 1999), but 
softer than the non-thermal power law component index 
derived by Choudhury et al.(2002) with ASM+BATSE 
data. 

4.3.3. EXO2030+375 




100 
Energy (keV) 



Fig. 10. The average spectrum of EXO 2030+375 for a 
total exposition time of 345 ks. The best fit power law 
(solid line) has an index of 2.72 ± 0.02 and the best fit 
thermal model (dotted line) has a temperature of 27.2 ± 
0.1. Upper limits are at 1 a. 

EXO 2030+375 is also detected up to 150 keV (fig. 
10). Its average photon spectrum (revolution 19 to 22, 
345 ks exposure) is well fitted by a thermal model with 
temperature ~ 27 keV or a power law of index ~ 2.7. 

A possible spectral feature at ~ 36 keV, tentatively ex- 
plained as a cyclotron absorption line, has been suggested 
by Reig and Coe (1999). 

We do not detect this feature in the averaged spectrum 
with a 2 a upper limit between 34 and 38 keV of 8 x 10 -5 
ph-cm -2 - s _1 . We cannot exclude the presence of such 
a feature in shorter timescales, but this requires further 
analysis. 

5. Discussion and Conclusion 

The emission from Cyg X-l has been detected up to 1 
MeV. While, in the energy range between 20 and 300 
keV, its flux varies by a factor of ~ 2, its spectral shape 
remains remarkably stable, on a one day timescale. Our 
mean photon spectrum is compatible with the average 
CGRO/OSSE spectrum in the hard state (McConnell et 
al., 2002). 

An outburst of EXO 2030+375 occured during our 
observations. The average SPI spectrum during this 6-7 
days outburst is compatible with that obtained by BATSE 
(Wilson et al, 2002). 



Initial results obtained on Cyg X-l, Cyg X-3 and EXO 
2030+375 are coherent with the knowledge we have on 
those sources and allow us to proceed with a more detailled 
analysis. 

Note that, although the angular distance between Cyg 
X-3 and EXO 2030+375 is close to the angular resolution 
of SPI, we can easily measure their fluxes and positions. 

Cyg X-l is one of the strongest persistent source in 
the SPI/INTEGRAL energy domain. From a technical 
point of view, it is used to test the coherence of our data 
analysis system, performance verification and inflight cal- 
ibrations. Those objectives have been achieved and they 
demonstrate that SPI is working perfectly. First of all, at 
the present time, the quality of the images leads us to con- 
clude that the spatial and energy responses are accurate 
to about 2-3 %. 

While qualitative results will not change, quantitative 
results are subject to the incertainties in the used energy 
response. We used the May 2003 version of the response 
matrix (Sturner et al., this issue), which underestimates 
the efficiency at low energy, resulting in steeper spectra. 

We plan to improve the data analysis system by using 
more precise background analysis and modelling and by 
including source variability in the image analysis package. 
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